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1 Introduction 

Particle correlations play an important role in the un¬ 
derstanding of multiparticle production. Correlations be¬ 
tween identical bosons, called Bose-Einstein correlations 
(BEC), are a well-known phenomenon in high-energy and 
nuclear physics (for reviews see PHI)- The BEC are of¬ 
ten considered to be the analogue of the Hanbury-Brown 
and Twiss effect in astronomy, describing the in¬ 

terference of incoherently emitted identical bosons ITfillTOl . 
They represent a sensitive probe of the space-time geom¬ 
etry of the hadronization region and allow the determina¬ 
tion of the size and the shape of the source from which 
particles are emitted. 

The production of identical bosons that are close to¬ 
gether in phase space is enhanced by the presence of BEC. 
The first observation of BEC effects in identically charged 
pions produced in pp collisions was reported in Ref. ESI 
[2T1 . Since then, BEC have been studied for systems of 
two or more identical bosons produced in various types of 
collisions, from leptonic to hadronic and nuclear collisions 
(see Refs. m and references therein). 

Studies of the dependence of BEC on particle mul¬ 
tiplicity and transverse momentum are of special inter¬ 
est. They help to understand the multiparticle production 
mechanism. The size of the source emitting the correlated 
particles has been observed to increase with particle multi¬ 
plicity. This can be understood as arising from the increase 
in the initial geometrical region of overlap of the colliding 
objects [22]: a large overlap implies a large multiplicity. 


While this dependence is natural in nucleus-nucleus colli¬ 
sions, the increase of size with multiplicity has also been 
observed in hadronic and leptonic interactions. In the lat¬ 
ter, it is understood as a result of superposition of many 
sources EE1EM271 or related to the number of jets 1281125] . 
High-multiplicity data in proton-proton interactions can 
serve as a reference for studies of nucleus-nucleus colli¬ 
sions. The effect is reproduced in both the hydrodynam- 
ical/hydrokinetic [30Tf32] and Pomeron-based [33JSJ ap¬ 
proaches for hadronic interactions where high multiplici¬ 
ties play a crucial role. The dependence on the transverse 
momentum of the emitter particle pair is another impor¬ 
tant feature of the BEC effect ES3. In nucleus-nucleus col¬ 
lisions the dependence of the particle emitter size on the 
transverse momentum is explained as a “collective flow”, 
which generates a characteristic fall-off of the emitter size 
with increasing transverse momentum |.'Hil l.'LSI while strong 
space-time momentum-energy correlations may offer an 
explanation in more “elementary” leptonic and hadronic 
systems [61171191 13011321135] where BEC measurements serve 
as a test of different models [5011521155H35] . 

In the present analysis, studies of one-dimensional BEC 
effects in pp collisions at centre-of-mass energies of 0.9 
and 7 TeV, using the ATLAS detector [37] at the Large 
Hadron Collider (LHC), are presented. At the LHC, BEC 
have been studied by the CMS [481149] and ALICE [501151] 
experiments. In the analysis reported here, the studies 
are extended to the region of high-multiplicities available 
thanks to the high multiplicity track trigger. The results 
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are compared to measurements at the same or lower en¬ 
ergies. 


2 Analysis 

2.1 Two-particle correlation function 


Bose-Einstein correlations are measured in terms of a two- 
particle correlation function, 


C 2 {Pl,P2) 


P(Pl,P2) 
Po(pi,Pz) ’ 


(1) 


where p\ and P2 are the four-momenta of two identical 
bosons in the event, p is the two-particle density function, 
and po is a two-particle density function (known as the 
reference function) specially constructed to exclude BEC 
effects. The densities p and po are normalized to unity, i.e. 
they are the probability density functions. 

In order to compare with data over the widest possi¬ 
ble range of centre-of-mass energies and system sizes, the 
density function is parameterized in terms of the Lorentz- 
invariant four-momentum difference squared, Q 2 , of the 
two particles, 

Q 2 = ~(Pi -P 2) 2 ■ (2) 

The BEC effect is usually described by a function with 
two parameters: the effective radius parameter R and the 
strength parameter A [55], where the latter is also called 
the incoherence or chaoticity parameter. A typical func¬ 
tional form is 

c 2(Q) = j^ = c o[i + n(\,QRm + £Q). (3) 


In a simplified scheme for fully coherent emission of identi¬ 
cal bosons, A = 0, while for incoherent (chaotic) emission, 
A = 1. The QR dependence comes from the Fourier trans¬ 
form of the distribution of the space-time points of bo¬ 
son emission. Several different functional forms have been 
proposed for L?(A, QR). Those used in this paper are de¬ 
scribed in Sect. m The fitted parameter e takes into ac¬ 
count long-distance correlations not fully removed from 
po- Finally, Cq is a normalization constant, typically cho¬ 
sen such that (72 (Q) is unity for large Q. In this paper, 
the density function p is calculated for like-sign charged- 

particle pairs, with both the ++ and-combinations 

included, p(Q) = p(++, -). All particles are treated as 

charged pions and no particle identification is attempted. 
The purity of the analysis sample in terms of identical bo¬ 
son pairs is estimated from MC to be about 70% (where 
about 69% are 7r ± 7r ± and about 1% are K^K^). The ef¬ 
fect of the purity is absorbed in the strength parameter 
A, while the results of the analysis on the effective radius 
parameter R were found to be not affected. 


2.2 Coulomb correction 

The long-range Coulomb force causes a momentum shift 
between the like-sign and unlike-sign pairs of particles. 


The density distributions are corrected for this effect by 
applying the Gamow penetration factor per track pair 
with a weight 1/G(Q) [55H55] (for review see Ref. [55]) 

PcoAQ) = , (4) 

where the Gamow factor G(Q) is given by 


G(Q) 


27r£ 

e 27rf _ \ 


(5) 


with the dimensionless parameter £ defined as 


am 

C= ~Q' 


( 6 ) 


Here a is the electromagnetic fine-structure constant and 
m is the pion mass. The sign of £ is positive for like- 
sign pairs and negative for unlike-sign pairs. The resulting 
correction on p(Q) decreases with increasing Q and at Q = 
0.03 GeV it is about 20%. A systematic uncertainty on 
G(Q) is considered to cover effects like the extended size 
of the emission source and other effects, see discussion in 
Ref. [TOFT] . Neither the Coulomb interaction nor the BEC 
effect are present in the generation of MC event samples 
which are used in the analysis. The Coulomb correction is 
thus not applied to MC events. 


2.3 Reference sample 

A good choice of the reference sample is important to al¬ 
low the experimental detection of the BEC signal. Ideally, 
Po(Q) should include all momentum correlations except 
those arising from BEC. Thus, several different choices 
have been studied to construct an appropriate reference 
sample. 

Most of the proposed approaches use random pairing 
of particles, such as mixing particles from different events 
(the “mixed event” technique [56]), or choosing them from 
the same event but from opposite hemispheres or by ro¬ 
tating the transverse momentum vector of one of the par¬ 
ticles of the like-sign pair [9]. Although these mixing tech¬ 
niques reproduce the topology and some properties of the 
event under consideration and destroy BEC, they vio¬ 
late energy-momentum conservation. Moreover, there are 
many possible ways to construct the pairs, such as mixing 
the particles randomly, or keeping some topological con¬ 
straints such as the event multiplicity, the invariant mass 
of the pair or the rapidity of the pair. All of these introduce 
additional biases in the BEC observables. For example, 
it was observed in dedicated MC studies that the single¬ 
ratio correlation functions C 2 using reference samples con¬ 
structed with the event mixing or opposite hemispheres 
techniques exhibit an increase in the low-Q BEC sensi¬ 
tive region. This effect is found to be more pronounced 
with increase of the multiplicity or average particle-pair 
transverse momentum and indicates that these reference 
samples are not suitable. 
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A natural choice is to use the unlike-sign particle pairs 
from the same events that are used to form pairs of like- 
sign particles, i.e., po(Q) = p(-\ —), called in the follow¬ 
ing the unlike-charge reference sample. This sample has 
the same topology and global properties as the like sign 
sample p(++, -), but is naturally free of any BEC ef¬ 

fect. Studying the C 2 correlation functions on MC, none of 
the deficits of the event mixing and opposite hemispheres 
techniques described above were observed. However, this 
sample contains hadron pairs from the decay of resonances 
such as p, 77 , 77 ', w, (f) , K*, which are not present in the 
like-sign combinations. These contribute to the low-Q re¬ 
gion and can give a spurious BEC signature with a large 
effective radius of the source Ema¬ 
il! this paper, the unlike-charge reference sample is 
used. To account for the effects of resonances, the two- 
particle correlation function (^(Q) is corrected using Monte 
Carlo simulation without BEC effects via a double-ratio 
I ?2 (Q) defined as 


MQ) 


C*(Q) 

Cf c {Q) 


?(++,—) / p mc (++,—) 
p(+~) / p mc (+~) 


2.4 The parameterizations of BEC 

Various parameterizations of the fi( A, QR) function can 
be found in the literature, each assuming a different shape 
for the particle-emitting source. In the studies presented 
here, the data are analysed using the following parame¬ 
terizations: 

— the Goldhaber parameterization prm] of a static Gaus¬ 
sian source in the plane-wave approach, 

17 = A • exp ( -R 1 2 Q 2 ) , ( 8 ) 

which assumes a spherical shape with a radial Gaus¬ 
sian distribution of the emitter; 

— the exponential parameterization of a static source 

17 = A • exp ( —RQ ), (9) 

which assumes a radial Lorentzian distribution of the 
source. This parameterization provides a better de¬ 
scription of the data at small Q values, as discussed 
in E]. 

The first moment of the fi(QR) distribution corresponds 
to 1 /R for the exponential form and to 1/ (Rt/tt) for the 
Gaussian form. To compare the values of the radius pa¬ 
rameters obtained from the two functions, the R value of 
the Gaussian should be compared to R/^/tt of the expo¬ 
nential form. 


3 Experimental details 

3.1 The ATLAS detector 

The ATLAS detector 07] is a multi-purpose particle physics 
experiment operating at one of the beam interaction points 


of the LHC. The detector covers almost the whole solid 
angle around the collision point with layers of tracking de¬ 
tectors, calorimeters and muon chambers. It is designed to 
study a wide range of physics topics at LHC energies. For 
the measurements presented in this paper, the tracking 
devices and the trigger system are of particular impor¬ 
tance. 

The innermost part of the ATLAS detector is the in¬ 
ner detector (ID), which has full coverage in <f> and covers 
the pseudorapidity range |r/| < 2.50 It consists of a silicon 
pixel detector (Pixel), a silicon microstrip detector (SCT) 
and a transition radiation tracker (TRT). These detectors 
are immersed in a 2 T solenoidal magnetic field. The Pixel, 
SCT, and TRT detectors have typical position resolutions 
of 10, 17 and 130 pm for the r-cf> coordinate, respectively. 
In the case of the Pixel and SCT, the resolutions are 115 
and 580 pm, respectively, for the second measured coor¬ 
dinate. A track from a charged particle traversing the full 
radial extent of the ID would typically have three Pixel 
hits, eight or more SCT hits and more than 30 TRT hits. 

The ATLAS detector has a three-level trigger system: 
Level 1 (LI), Level 2 (L2) and Event Filter (EF). For this 
measurement, the trigger relies on the LI signals from the 
Beam Pickup Timing devices (BPTX) and the Minimum- 
Bias (MB) Trigger Scintillators (MBTS). The BPTX are 
composed of electrostatic button pick-up detectors attached 
to the beam pipe and located 175 m from the centre of the 
ATLAS detector in both directions along the beam pipe. 
The MBTS are mounted at each end of the detector in 
front of the liquid-argon end-cap calorimeter cryostats at 
z = ±3.56 m. They are segmented into eight sectors in az¬ 
imuth and two rings in pseudorapidity (2.09 < \q\ < 2.82 
and 2.82 < \r]\ < 3.84). Data was collected requiring co¬ 
incidence of BPTX and MBTS signals, where only a single 
hit in the MBTS was required on either side of the detec¬ 
tor. The efficiency of this trigger was studied with events 
collected with a separate prescaled LI BPTX trigger, fil¬ 
tered by ID requirements at L2 and at EF level in order 
to obtain inelastic interactions and found to be 98% for 
two selected tracks and 100 % for more than four selected 

tracks [MIES- 

High-multiplicity track (HM) events were collected at 
7 TeV using a dedicated high-multiplicity track trigger. At 
LI, the collisions were triggered using the summed trans¬ 
verse energy (EEt) in all calorimeters, calibrated at the 
electromagnetic energy scale | 66 j. The high-multiplicity 
events were required to have EEt > 20 GeV. A high num¬ 
ber of hits in the SCT was required at L2, while at the EF 
level at least 124 tracks with px > 400 MeV were required 
to originate from a single vertex. 


1 ATLAS uses a right-handed coordinate system with its ori¬ 

gin at the nominal interaction point (IP) in the centre of the 
detector and the 2 -axis along the beam pipe. The i-axis points 
from the IP to the centre of the LHC ring, and the j/-axis 
points upward. Cylindrical coordinates (r, tj >) are used in the 
transverse plane, <j> being the azimuthal angle around the beam 
pipe. The pseudorapidity is defined in terms of the polar angle 
(9 as p = — lntan(0/2). 
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3.2 Data and Monte Carlo samples 

The study is carried out using the pp-collision datasets at 
the centre-of-mass energies y/s = 0.9 and 7 TeV that were 
used in previously published ATLAS studies of minimum- 
bias interactions [MM5] . 

The event and track selection criteria are the same 
as the ones used for the ATLAS minimum-bias multiplic¬ 
ity analysis [B5] with the same minimum-bias trigger and 
quality criteria for the track reconstruction. All events in 
these datasets are required to have at least one vertex mi 
formed from a minimum of two tracks with px > 100 MeV 
and consistent with the average beam spot position within 
the ATLAS detector (primary vertex) |B5]. The tracks sat¬ 
isfying the above-mentioned selection criteria are used as 
the input to determine the corrected distributions, as de¬ 
scribed in Sect. The multiplicity of selected tracks 
with px > 100 MeV and |p| < 2.5 within an event is de¬ 
noted by n se \. 

The contributions from beam-gas collision and from 
non-collision background (cosmic rays and detector noise) 
were investigated in Ref. [OH and found to be negligible. 
Events with more than one primary vertex (less than 0.3% 
of the sample) are rejected in order to prevent a bias from 
multiple proton-proton interactions (pile-up) in the col¬ 
liding proton bunches. 

The same event selection criteria are applied to high- 
nrultiplicity events, which are defined to be those with 
at least 120 selected tracks. To estimate the possible in¬ 
fluence of multiple pp interactions in the 7 TeV high- 
nrultiplicity track trigger data, the distribution of the dis¬ 
tances Az between the z coordinates of primary and pile- 
up vertices are studied. The study shows that on average 
there is less than one pile-up track selected in the HM 
sample, which has a negligible influence on the BEC stud¬ 
ies. 

For the measurements at y/s = 0.9 TeV, about 3.6 x 10 5 
events with a total of more than 4.5 x 10 6 tracks are after 
selection, and in the case of /s = 7 TeV, about 10' events 
with about 2.1 x 10 8 tracks overall are after selection. This 
corresponds to integrated luminosities of ~7 pb -1 and 
~190 pb _1 at 0.9 TeV and 7 TeV, respectively. For the 
measurements at 7 TeV with the high-nrultiplicity track 
trigger, about 1.8 x 10 4 events with more than 2.7 x 10 6 
tracks overall were after selection. This corresponds to in¬ 
tegrated luminosity of ~12.4 nb _1 . 

Large Monte Carlo samples of minimum-bias and high- 
multiplicity events were generated using the PYTHIA 6.421 
Monte Carlo event generator [69] with the ATLAS MC09 
set of optimised parameters (tune) [7U] ( 1-1 x 10 7 for y/s = 
900 GeV, 2.7 x 10 7 for /I = 7 TeV and 1.8 x 10 6 for 
/s = 7 TeV high-multiplicity data) with non-diffractive, 
single-diffractive and double-diffractive processes included 
in proportion to the cross sections predicted by the model. 
As discussed in Sec. 12.21 no simulation of the BEC ef¬ 
fect is implemented in the generator. This is the baseline 
Monte Carlo generator which reproduces single-particle 
spectra [Ml- The generated events were passed through 
the ATLAS simulation and reconstruction chain; the de¬ 
tector simulation program [71] is based on GEANT4 [72 . 


Dedicated sets of high-nrultiplicity events were also gen¬ 
erated. 

For the study of systematic effects, additional Monte 
Carlo samples were produced using the PHOJET 1.12.1.35 
generator [73] , PYTHIA with the PerugiaO tune [71] ; and 
the EPOS 1.99_v2965 generator [TS] for the high-multiplici¬ 
ty analysis. The PHOJET program uses the Dual Par- 
ton Model m for low-px physics and is interfaced to 
PYTHIA for the fragmentation of partons. The EPOS 
generator is based on an implementation of the QCD- 
inspired Gribov-Regge field theory describing soft and 
hard scattering simultaneously, and relies on the same par- 
ton distribution functions as used in PYTHIA. The EPOS 
LHC tune is used with parameters optimised to describe 
the LHC minimum-bias data I76j . 

The high-nrultiplicity PYTHIA MC09 and EPOS sam¬ 
ples, each are about two magnitudes larger than the data 
sample. The C 2 (Q) single-ratio correlation functions in 
MC reproduce data well for Q > 0.5 GeV. In the region 
Q < 0.5 GeV, the BEC effect is clearly seen in the data 
C 2 (Q) correlation function while no such effect is seen in 
the MC as expected, since no BEC present in MC. 


3.3 Data correction procedure 

Following the procedure applied in the previous ATLAS 
minimum-bias measurements I7RTK51 , each track is assigned 
a weight which corrects for the track reconstruction effi¬ 
ciency, for the fraction of secondary particles, for the frac¬ 
tion of the primary particles^ outside the kinematic range 
and for the fraction of fake tracks0 In addition, the effect 
of events lost due to trigger and vertex reconstruction in¬ 
efficiencies is corrected for using an event-by-event weight 
applied to pairs of particles in the Q distribution. The effi¬ 
ciency of the high-nrultiplicity track trigger has been stud¬ 
ied in data as a function of the number of reconstructed 
tracks and is found to be 5% for 120 selected tracks and 
to reach a plateau at 100% once 150 tracks are selected. 
The measured trigger inefficiency is used to correct the 
experimental distributions and is found to have negligible 
impact on the extraction of the BEC parameters discussed 
in Sec. [5] 

The multiplicity distributions are corrected to the par¬ 
ticle level using an iterative method that follows the Bayesi¬ 
an approach m as it is described in Refs. 1641165] . An un¬ 
folding matrix reflecting the probability of reconstructing 
n S ei charged tracks in an event with generated charged- 
particle multiplicity n c h is populated using Monte Carlo 
simulation and applied to the data. The unfolding matrix 
is built using the ATLAS MC09 PYTHIA tune [TU]. The 
unfolding procedure converges after the fifth iteration. It 
is found that the corrected multiplicity distribution agrees 

2 In the Monte Carlo simulations, primary charged particles 
are defined as charged particles with a mean lifetime r > 0.3 x 
10~ 10 s either directly produced in pp collisions or from the 
subsequent decay of particles with a shorter lifetime. 

3 Fake tracks are tracks constructed from tracker noise 
and/or hits which are not produced by a single-particle. 
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well with the published result [641165] . The unfolding pro¬ 
cedure of the 7 TeV high-multiplicity data follows the same 
technique and unfolding matrix used in the previous anal¬ 
ysis of minimum-bias data in Ref. [M], restricted to the 
region of high charged particle multiplicity specific to this 
analysis, and convolved with a normalised Gaussian distri¬ 
bution to account for the experimental resolution on the 
number of selected tracks. It is found that a number of 
120 selected tracks at detector level, n se i, corresponds to 
about 150 charged particle, n c h, at particle level. Momen¬ 
tum distributions are unfolded in a similar way. 

For all distributions, closure tests are carried out us¬ 
ing Monte Carlo samples corrected according to the same 
procedure as used in the data. The difference obtained be¬ 
tween the reweighted distributions and those at the par¬ 
ticle level is due to tracking effects such as a smaller re¬ 
construction efficiency for pairs of tracks with very small 
opening angle. These effects are small for correlation func¬ 
tions constructed using data, typically 1-3%, and are in¬ 
cluded in the systematic uncertainty. In the case of the 
unfolded Q distributions, the data are corrected for the 
bias from secondary tracks using Monte Carlo simulation 
and the corresponding systematic uncertainty is obtained 
by variation of the amount of material in the inner detec¬ 
tor by ±10%. 

4 Systematic uncertainties 

The systematic uncertainties of the inclusive fit parame¬ 
ters, R and A, of the exponential model are summarized 
in Table [Q The following contributions to the systematic 
uncertainties on the fitted parameters are considered. 

The systematic uncertainties resulting from the track 
reconstruction efficiency, which are parameterized in bins 
of px and 77 , were determined in earlier analyses [641165] . 
These cause uncertainties in the track weights of parti¬ 
cle pairs in the Q distributions entering the correlation 
functions. 

The effects of track splitting and merging are sizeable 
only for very low Q values (smaller than 5 MeV), and are 
found to be negligible for the measurements with Q > 20 
MeV. 

The leading source of systematic uncertainty is due to 
differences in the Monte Carlo generators used to calcu¬ 
late the i ?2 correlation function from the C 2 correlation 
function. The corresponding contribution to the system¬ 
atic uncertainty is estimated as the root-mean-squared 
(RMS) spread of the results obtained for the different 
Monte Carlo datasets. The statistical uncertainties aris¬ 
ing from the Monte Carlo datasets are negligibly small. 

The systematic uncertainty due to Coulomb correc¬ 
tions is estimated by varying the corrections by ±20%. 

The influence of the fit range is estimated by changing 
the upper bound of the Q range from the nominal 2 GeV: 
decreasing it to 1.5 GeV and increasing it up to 2.5 GeV. 
The latter better estimates the uncertainty due the long- 
range correlations. This contribution is taken into account 
by the value of e, the parameter in the linear term of 
Eq. ([3]) describing the long-range correlations. 


Other effects contributing to the systematic uncertain¬ 
ties are the lowest value of Q for the fit, the bin size and 
exclusion of the interval 0.5 < Q < 0.9 GeV due to the 
overestimate of the p meson contribution in the Monte 
Carlo simulations, as discussed in the following Sec. 15.11 
These uncertainties are estimated by varying the lowest 
Q value in the fit by ±10 MeV, by changing the bin size 
by ±10 MeV, and by broadening the excluded interval by 
100 MeV on both sides. 

The background of photon conversions into e + e - pairs 
was studied and found to be negligible. 

To test the effect of treating all charged particles as 
pions, the double-ratio correlation functions i ?2 are also 
obtained using only identical particles in the Monte Carlo 
sample to compute the correction. The resulting BEC pa¬ 
rameters fitted to the i ?2 functions defined this way show 
negligible differences to the nominal result and no further 
systematic uncertainties are assigned. 

Finally, the systematic uncertainties are combined by 
adding them in quadrature and the resulting values are 
given in the bottom row of Table [T] 

The same sources of uncertainty are considered for the 
differential measurements in n c h and the average trans¬ 
verse momentum fcx of a pair, and their impact on the fit 
parameters is found to be similar in size. 
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Table 1 . Systematic uncertainties on A and R for the exponential fit of the two-particle double-ratio correlation function 
Ri{Q) in the full kinematic region at y/s =0.9 and 7 TeV for minimum-bias and high-multiplicity events. 



0.9 TeV 

7 TeV 

7 TeV (HM) 

Source 

A 

R 

A 

R 

A 

R 

Track reconstruction efficiency 

0 .6% 

0.7% 

0.3% 

0 .2% 

1.3% 

0.3% 

Track splitting and merging 

negligible 

negligible 

negligible 

Monte Carlo samples 

14.5% 

12.9% 

7.6% 

10.4% 

5.1% 

8.4% 

Coulomb correction 

2 .6% 

0 .1% 

5.5% 

0 .1% 

3.7% 

0.5% 

Fitted range of Q 

1 .0% 

1 .6% 

1 .6% 

2 .2% 

5.5% 

6 .0% 

Starting value of Q 

0.4% 

0.3% 

0.9% 

0 .6% 

0.5% 

0.3% 

Bin size 

0 .2% 

0 .2% 

0.9% 

0.5% 

4.1% 

3.4% 

Exclusion interval 

0 .2% 

0 .2% 

1 % 

0 .6% 

0.7% 

1 .1% 

Total 

14.8% 

13.0% 

9.6% 

10.7% 

9.4% 

10.9% 


5 Results 

5.1 Two-particle correlations 

In Fig. the double-ratio R 2 (Q) distributions, measured 
for 0.9 and 7 TeV, are compared with Gaussian and ex¬ 
ponential fitting functions, Eqs. ([5]) and ([9]). The fits are 
performed in the Q range 0.02 GeV to 2 GeV and with a 
bin width of 0.02 GeV. The upper Q limit is chosen to be 
far away from the low-Q region, which is sensitive to BEC 
effects and resonances. Around Q ~ 0.7 GeV there is a vis¬ 
ible bump which is due to an overestimate of p — > 7r + 7r - 
decays in the Monte Carlo simulation. Therefore the re¬ 
gion 0.5 < Q < 0.9 GeV is excluded from the fits. As seen 
in Fig. [I] the Gaussian function does not describe the low- 
13 region while the exponential function provides a good 
description of the data. 

The resolution of the Q variable is better than 10 MeV 
for the region most sensitive to BEC effect, Q < 0.4 GeV. 
The Q resolution is included in the fit of i? 2 by convolving 
the fitting function with a Gaussian detector resolution 
function. The change in the fit results from those with no 
convolution applied is found to be negligible. 

In the process of fitting i? 2 (Q) with the exponential 
function, large y 2 values are observed, in particular for 
the 7 TeV sample where statistical uncertainties on the 
fitted data points are below 2-4%. These large y 2 val¬ 
ues can be traced back to a small number of individual 
points or small cluster of points. The removal of these 
points does not change the results of the fit while the y 2 
substantially improves. In the analysis of the 7 TeV data, 
for most of the considered cases, the expected statistical 
uncertainties are small compared to the systematic ones, 
therefore only total uncertainties on the fitted parameters 
are given. The latter include the statistical uncertainties 
rescaled by ^/y 2 /ndf 175]. For consistency, the same treat¬ 
ment is applied to the 0.9 TeV analysis where the statis¬ 
tical uncertainties are of the same order of magnitude as 
the systematic ones. 

The results of BEC parameters for exponential fits of 
the two-particle double-ratio correlation function i? 2 (Q) 
for events with the unlike-charge reference sample are 

A = 0.74 ± 0.11, R = (1.83 ± 0.25) fm at y/s = 0.9 TeV for n ch > 2, 


A = 0.71 ± 0.07, R = (2.06 ± 0.22) fmat\/s = 7 TeV for n c h > 2, 

A = 0.52 ± 0.06, R = (2.36 ± 0.30) fm at y/s = 7 TeV for n ch > 150. 

The values of the fitted parameters are close to the 
values obtained by the CMS [49] and ALICE [50] experi¬ 
ments. 


5.2 Multiplicity dependence 

The i? 2 (Q) functions defined in Eq. J7]), are shown for var¬ 
ious multiplicity intervals in Fig.[2]for 0.9 TeV, 7 TeV and 
7 TeV high-multiplicity data. The multiplicity intervals 
are chosen so as to be similarly populated and compara¬ 
ble to those used by other LHC experiments [48U5T] . Only 
the exponential fit is shown. As in the fit procedure for the 
inclusive case, the detector Q resolution is included in the 
fits. 

Within the multiplicity studies, the BEC parameters 
are also measured by excluding the low-multiplicity events, 
7i c h < 8, expected to be contaminated by diffractive physics 
m- No noticeable changes in the strength and radius pa¬ 
rameters for n c h > 8 are observed compared to the full 
multiplicity range for n c h > 2. 

The multiplicity dependence of the A and R parame¬ 
ters is shown in Fig. [3] The A parameter decreases with 
multiplicity, faster for 0.9 TeV than for 7 TeV interac¬ 
tions. The decrease of the A parameter with n c h is found 
to be well fitted with the exponential function A(n c h) = 
,y e -<5nch_ The fit parameter values are presented in Table 
[2] for 0.9 TeV and for the combined nominal and high- 
multiplicity 7 TeV data. 

The R parameter increases with multiplicity up to 
about rich — 50 independently of the center of mass en¬ 
ergy. For higher multiplicities, the measured R parameter 
is observed to be independent of multiplicity. For n c h < 82 
at 0.9 TeV and n c h < 55 at 7 TeV the n c h dependence of 
R is fitted with the function R(n c h) = ot^/n c h, similar to 
that used in heavy-ion studies EHEP- The results of the 
fit are presented in Table [2] and are close to the CMS re¬ 
sults [32]. The fit parameters do not change significantly 
within uncertainties if data points with n c h > 55 are in¬ 
cluded in the fit, while the quality of the fit significantly 
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Fig. 1. The two-particle double-ratio correlation function R^iQ) for charged particles in pp collisions at (a) y/s =0.9 TeV, (b) 
7 TeV and (c) 7 TeV high-multiplicity events. The lines show the Gaussian and exponential fits as described in the legend. The 
region excluded from the fits is indicated. The error bars represent the statistical uncertainties. 


degrades. Therefore the fit is limited to the data points 
with n c h < 55. The n c h dependence of R at 7 TeV is fit¬ 
ted with a constant R(n c h) = /? for n c h > 55; the resulting 
value is given in Table [2] Qualitatively CMS [49] and UAl 
El results for the radius parameter follow the same trend 
as a function of n c h as ATLAS data points up to n c h < 55. 
The ATLAS and ALICE [5(711511 results on the multiplic¬ 
ity dependence of the radius parameter cannot be directly 
compared due to much narrower rj region used by ALICE. 

The observed change of the fitted parameters with 
multiplicity has been predicted in Refs. [51 I25H2?] . and is 
similar to the one also observed in e + e _ interactions [251 . 
however the saturation of R for very high multiplicity is 
observed for the first time. 

The saturation of R at high multiplicities is expected in 
a Pomeron-based model PIIM] as the consequence of the 
overlap of colliding protons, with the value of the radius 


parameter at n c h ~ 70 close to the one obtained in the 
present studies. However, the same model predicts that 
above n c h « 70, R will decrease with multiplicity, return¬ 
ing to its low-multiplicity value which is not supported by 
the data. 


5.3 Dependence on the transverse momentum of the 
particle pair 

The average transverse momentum fcx of a particle pair is 
defined as half of the magnitude of the vector sum of the 
two transverse momenta, fcx = |Pt i +Pt,2|/2. The study 
is performed in the fcx intervals which are chosen in a 
way to be similarly populated and, as for the multiplicity 
bins, to be similar to the intervals used by other LHC 
experiments I48H5T] . 
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(b) 



(c) 

Fig. 2. The two-particle double-ratio correlation function R- 2 (Q) for charged particles in pp collisions for multiplicity intervals 
(a) 36 < n c h < 45 at y/s =0.9 TeV, (b) 68 < n c h < 79 at 7 TeV and (c) 183 < n c h < 197 at 7 TeV high-multiplicity events. 
The lines show the results of the exponential fit. The region excluded from the fits is indicated. The error bars represent the 
statistical uncertainties. 


Table 2. Results of fitting the multiplicity, n c h, and the transverse momentum of the pair, Ict, dependence of the BEC 
parameters R and A with different functional forms and for different data samples. The error represent the quadratic sum of 
the statistical and systematic uncertainties. 


BEC 
par am. 

Fit 

function 


0.9 TeV 


7 TeV 


Minimum-bias events 

High-multiplicity events 

-^(^ch) 

a-^n^h 

a 

= 0.64 ± 0.07 fm (n ch < 82) 

a = 

0.63 ± 0.05 fm (n ch < 55) 


P 


— 


P = 2.28 ± 0.32 fm (n ch > 55) 

A (rich) 

cyg-^ch 

7 

= 1.06 ±0.10 


7 = 0.96 ±0.07 



S 

= 0.011 ±0.004 


S = 0.0038 ± 0.0008 

R{k T ) 

£ e ~Kk T 

e 

= 2.64 ± 0.33 fm 

f = 

2.88 ± 0.27 fm 

£ = 3.39 ± 0.54 fm 



K, 

= 1.48 ± 0.67 GeV” 1 

K = 

1.05 ± 0.58 GeV -1 

k = 0.92 ± 0.73 GeV -1 

A(fcr) 

pe~^ 


= 1.20 ±0.18 

n = 

1.12 ±0.10 

p = 0.75 ±0.10 



V 

= 2.00 ± 0.35 GeV" 1 

V = 

1.54 ±0.26 GeV -1 

v = 0.91 ± 0.45 GeV' 1 
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Fig. 3. Multiplicity, n c h, dependence of the parameters (a) A and (b) R obtained from the exponential fit to the two-particle 
double-ratio correlation functions R 2 (Q) at yfs =0.9 and 7 TeV. The solid and dashed curves are the results of (a) the exponential 
and (b) ^/n c h for n c h < 55 fits. The dotted line in (b) is a result of a constant fit to minimum-bias and high-multiplicity events 
data at 7 TeV for n c h > 55. The error bars represent the quadratic sum of the statistical and systematic uncertainties. 
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As an example, the i? 2 (Q) distributions for the 500 < 
fcx < 600 MeV interval for the 0.9 TeV, 7 TeV and high- 
multiplicity 7 TeV samples are shown in Fig. [1] together 
with the results of the corresponding exponential fit. For 
the i? 2 (Q) correlation function measured at 7 TeV (see 
Fig. mb)), there is an indication that the Monte Carlo 
simulation overestimates the production and decay of the 
w-meson in the Q region of 0.3-0.44 GeV. This region is 
thus excluded from the fit range for fcx > 500 MeV bin 
results. 

In the region most important for the BEC parameters, 
the quality of the exponential fit is found to deteriorate as 
fcx increases. This is due to the fact that at large kf val¬ 
ues, the characteristic BEC peak becomes steeper than the 
exponential function can accommodate. Despite the dete¬ 
riorating fit quality, the behaviour of the fitted parameters 
is presented for comparison with previous experiments. 

The fit values of the A and R parameters are shown 
in Fig. [5] as a function of /ct- The values of both A and R 
decrease with increasing k t- 

The decrease of A with is well described by an expo¬ 
nential function, A(fcx) = ^e -J ' fcT . The fcx dependence of 
the R parameter is also found to follow an exponential de¬ 
crease, i?(fcx) = £ e~ KkT . The shapes of the fcx dependence 
are similar for the 7 TeV and the 7 TeV high-multiplicity 
data. The results of the fits are presented in Table [2j 

In Fig. [5](b), the fcx dependence of the R parameter 
is compared to the measurements performed by the E735 
[SU] and the STAR [ST] experiments with mixed-event ref¬ 
erence samples. These earlier results were obtained from 
Gaussian fits to the single-ratio correlation functions and 
therefore the values of the measured radius parameters 
are multiplied by as discussed in Sect. 12.41 The values 
of the parameters are observed to be energy-independent 
within the uncertainties. 

In Fig. O the fcx dependence of A and R, obtained 
for the 7 TeV data, is also studied in various multiplicity 
regions: 2 < n c h < 9; 10 < n c h < 24; 25 < n c h < 80; 
and 81 < n c h < 125. The decrease of A with fcx is nearly 
independent of multiplicity for n c h > 9 and the same as 
for the inclusive case. For n c h < 9 no conclusions can be 
drawn due to the large uncertainties. The i?-parameter 
decreases with fcx and exhibits an increase with increasing 
multiplicity as was observed for the fully inclusive case. 


The ATLAS Collaboration: Two-particle Bose-Einstein correlations 


11 





Fig. 4. The two-particle double-ratio correlation function Ri{Q ) for charged particles in pp collisions for 500 < < 600 

MeV interval at (a) y/s =0.9 TeV, (b) 7 TeV and (c) 7 TeV high-multiplicity events. The average transverse momentum Ict of 
the particle pairs is defined as hr = |pt,i + Pt,2|/2. The lines show the exponential fits. The region excluded from the fits is 
indicated. The error bars represent the statistical uncertainties. 
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Fig. 5. The k t dependence of the fitted parameters (a) A and (b) R obtained from the exponential fit to two-particle double¬ 
ratio at yfs =0.9 TeV, 7 TeV and 7 TeV high-multiplicity events. The average transverse momentum of the particle pairs 
is defined as /ct = |pt,i + Pt,2|/2. The solid, dashed and dash-dotted curves are results of the exponential fits for 0.9 TeV, 
7 TeV and 7 TeV high-multiplicity data, respectively. The results are compared to the corresponding measurements by the E735 
experiment at the Tevatron ESI . and by the STAR experiment at RHIC EH- The error bars represent the quadratic sum of the 
statistical and systematic uncertainties. 




Fig. 6. The /ct dependence of the fitted parameters (a) A and (b) R obtained from the exponential fit to the two-particle 
double-ratio correlation function i? 2 (Q) at yfs = 7 TeV for the different multiplicity regions: 2 < n c h < 9 (circles), 10 < n c h < 24 
(squares), 25 < n c h < 80 (triangles) and 81 < n c h < 125 (inverted triangles). The average transverse momentum of the 
particle pairs is defined as = |pt,i + Pt,2|/2. The error bars represent the quadratic sum of the statistical and systematic 
uncertainties. 
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6 Summary and conclusions Acknowledgments 


The two-particle Bose-Einstein correlations of like-sign 
hadrons with px > 100 MeV and \r/\ < 2.5 produced 
in pp collisions recorded by the ATLAS detector at 0.9 
and 7 TeV at the CERN Large Hadron Collider are stud¬ 
ied. In addition to minimum-bias data, high-multiplicity 
data recorded at 7 TeV using a dedicated trigger are in¬ 
vestigated. The integrated luminosities are about 7 /ib _1 , 
190 p b -1 and 12.4 nb -1 for 0.9 TeV, 7 TeV minimum-bias 
and 7 TeV high-multiplicity data samples, respectively. 

The studies were performed using the double-ratio cor¬ 
relation function. In the double-ratio method, the single¬ 
ratio correlation function obtained from the data is di¬ 
vided by a similar single-ratio calculated using Monte Carlo 
events, which do not have BEC effects. The reference sam¬ 
ple for each of the two single-ratios is constructed from 
unlike-sign charged-particle pairs. 

A clear signal of Bose-Einstein correlations is observed 
in the region of small four-momentum difference. To quan¬ 
titatively characterize the BEC effect, Gaussian and expo¬ 
nential parametrizations are fit to the measured correla¬ 
tion functions. As observed in studies performed by other 
experiments, the Gaussian parameterization provides a 
poor description of the BEC-enhanced region and hence 
the exponential parameterization is used for the final re¬ 
sults. 

The BEC parameters are studied as a function of the 
charged-particle multiplicity and the transverse momen¬ 
tum of the particle pair. A decrease of the correlation 
strength A along with an increase of the correlation source 
size parameter R are found with increasing charged-particle 
multiplicity. On the other hand no dependence of R on the 
centre-of-mass energy of pp collisions is observed. For the 
first time a saturation of the source size parameter is ob¬ 
served for multiplicity n c h > 55. The correlation strength 
A and the source size parameter R are found to decrease 
with increasing average transverse momentum of a pair. 
The study of BEC in (n c h, fcx) bins at 7 TeV shows a de¬ 
crease of the R parameter with k x for different multiplicity 
ranges, while the R values increase with multiplicity. The 
A parameter is found to decrease with fcx independently 
of the multiplicity range. These resemble the dependences 
for the inclusive case at 7 TeV for minimum-bias and high- 
multiplicity data. 

A comparison is made to the measurements by other 
experiments at the same and lower energies where pos¬ 
sible. The measurements presented here complement the 
earlier measurements by extending the studies to higher 
multiplicities and transverse momenta. This has allowed 
a first observation of a saturation in the magnitude of the 
source radius parameter at high charged-particle multi¬ 
plicities, and confirms the exponential decrease, observed 
in previous measurements of the radius parameters with 
increasing pair transverse momenta. 
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